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Abstract —In this paper, we address the problem of optimal 
power allocation at the relay in two-hop secure communications 
under practical conditions. To guarantee secure communication 
during the long-distance transmission, the massive MIMO (M- 
MIMO) relaying techniques are explored to significantly enhance 
wireless security. The focus of this paper is on the analysis and 
design of optimal power assignment for a decode-and-forward 
(DF) M-MIMO relay, so as to maximize the secrecy outage 
capacity and minimize the interception probahility, respectively. 
Our study reveals the condition for a nonnegative the secrecy 
outage capacity, obtains closed-form expressions for optimal 
power, and presents the asymptotic characteristics of secrecy 
performance. Finally, simulation results validate the effectiveness 
of the proposed schemes. 

I. Introduction 

The open nature of the wireless channel facilitates a mul¬ 
tiuser transmission, but also incurs the security problem. Re¬ 
cently, as a complement of traditional upper-layer encryption 
techniques, physical layer security (PHY-security) has been 
proposed to realize secure communications by making use of 
the characteristics of wireless channels, i.e., noise, fading and 
interference IT]. 

From an information-theoretic viewpoint, the performance 
of PHY-security is determined by the rate difference be¬ 
tween the legitimate channel and the eavesdropper channel 
m El. Therefore, to enhance wireless security, it makes 
sense to simultaneously increase the legitimate channel rate 
and decrease the eavesdropper channel rate. Inspired by this, 
various physical layer techniques have been introduced to 
improve the secrecy performance. Wherein, MIMO relaying 
technique gains considerable attention due to the following 
two reasons. First, the relay can provide a diversity gain 
and shorten the accessing distance, and thus improve the 
secrecy performance 0. Second, MIMO techniques, such as 
spatial beamforming, can reduce the information leakage to 
the eavesdropper 0. The beamforming schemes at the MIMO 
relay based on amplify-and-forward (AF) and decode-and- 
forward (DF) relaying protocols were presented in 0 and 
0, respectively. It is worth pointing out that the optimal beam 
design at the relay requires global channel state information 
(CSI) 0. Yet, the CSI, especially eavesdropper CSI is difficult 
to obtain, since the eavesdropper is usually passive and keeps 
silent. Therefore, it is impossible to realize absolutely secure 
communications over fading channels. In this context, statisti¬ 
cal secrecy performance metrics, e.g., secrecy outage capacity 
and interception probability, are adopted to evaluate wireless 


security 0. 

Recently, an advanced MIMO relaying technology, namely 
massive MIMO (M-MIMO) relaying, is introduced into secure 
communications to further improve the secrecy performance 
na. Even without eavesdropper CSI, M-MIMO relaying can 
generate a very high-resolution spatial beam, and thus the 
information leakage to the eavesdropper is quite small. More 
importantly, the secrecy performance can be enhanced by 
simply adding more antennas. Hence, the challenging issue 
of short-distance interception in secure communications can 
be well solved. Note that in two-hop secure communications, 
the transmit power at the relay has a great impact on the 
secrecy performance, since it affects the signal quality at 
both the destination and the eavesdropper. An optimal power 
allocation scheme for a multi-carrier two-hop single-antenna 
relaying network was given in by maximizing the sum 
secrecy rate. However, the power allocation for a multi-antenna 
relay, especially an M-MIMO relay, is still an open issue. In 
this paper, we focus on power allocation for DF M-MIMO 
secure relaying systems under very practical assumptions, 
i.e., no eavesdropper CSI and imperfect legitimate CSI. The 
contributions of this paper are three-fold; 

1) We reveal the relation between the secrecy outage capac¬ 
ity and a defined relative distance-dependent path loss, 
and then derive the condition for a nonnegative secrecy 
outage capacity and the constraint on the minimum 
number of antennas. 

2) We derive closed-form expressions for the optimal 
power at the relay in the sense of maximizing the 
secrecy outage capacity and minimizing the interception 
probability, respectively. 

3) We present clear insights into the secrecy performance 
through asymptotic analysis. We show that the maximum 
secrecy outage capacity is an increasing function of the 
source power while the minimum interception probabil¬ 
ity keeps constant. 

The rest of this paper is organized as follows. We first 
give an overview of the DF M-MIMO secure relaying system 
in Section II, and then derive two optimal power allocation 
schemes for the relay in Section III. In Section IV, we present 
simulation results to validate the effectiveness of the proposed 
schemes. Finally, we conclude the paper in Section V. 


II. System Model 


and 


We consider a time division duplex (TDD) two-hop massive 
MIMO (M-MIMO) relaying system, where a single antenna 
source transmits message to a single antenna destination with 
the aid of a relay with Nji antennas, while a passive single 
antenna eavesdropper intends to intercept the information. 
Note that the number of antennas at the relay in such an M- 
MIMO relaying system is quite large, i.e., Nji = 100 or even 
bigger. It is assumed that there is no direct transmission be¬ 
tween the source and the destination due to a long propagation 
distance. The relay system works in a half-duplex mode, which 
means any successful transmission requires two time slots. 
Specifically, the source sends the signal to the relay in the 
first time slot, and then the relay forwards the post-processing 
signal to the destination in the second time slot. In this paper, 
we assume the eavesdropper is far from the source and close 
to the relay, since it was concluded that the signal comes from 
the relay directly. Therefore, it is reasonable to assume that the 
eavesdropper only monitors the transmission from the relay to 
the destination. 

We use ,yas,Rhs,R, y/aR^D^R.D and y/aR^E^R,E to rep¬ 
resent the channels from the source to the relay, the relay to 
the destination, and the relay to the eavesdropper, respectively, 
where as,R, aR,D and ur^e are the distance-dependent path 
losses and hs,R, h/j d, and h/{ ^ denote the channel fading 
coefficient vectors with independent and identically distributed 
(i.i.d.) zero mean and unit variance complex Gaussian entries. 
It is assumed that the channels remain constant during a time 
slot and fade independently over slots. Thus, the received 
signal at the relay in the first time slot can be expressed as 

yR = V Psois,r^s,rs + n/j, (1) 

where s is the normalized Gaussian distributed transmit sym¬ 
bol, Ps is the transmit power at the source, and n^j stands 
for the additive white Gaussian noise with unit variance at the 
relay. We assume that the relay has perfect CSI about 
by channel estimation. Then, maximum ratio combination 
(MRC) decoding is performed to recover the information. 
Specifically, the received signal is multiplied by a vector 

During the second time slot, the relay forwards the decoded 
signal s through maximum ratio transmission (MRT). We 
assume that the relay has partial CSI about h/j D due to 
channel reciprocity impairment in TDD systems. The relation 
between the estimated CSI \iR E and the real CSI Iir^e is 
given by 

hfl.D = ^/f^R,D + \/l — pe, (2) 

where e is the error vector with i.i.d. zero mean and unit 
variance complex Gaussian entries, and is independent of 
hfl.D- P, scaling from 0 to 1, is the correlation coefficient 
between Iir^d and h/j £). Thus, the received signals at the 
destination and the eavesdropper are given by 

Vd = \J Pro-r^d^^^d^ + (3) 


UE = \J PROtR^E^^^E^ + ^E, (4) 

respectively, where Pr is the transmit power of the relay, r = 
\rs is the forwarded signal with \r = hfl,z)/||h_R,D|| being 
a MRT beamforming vector, no and are additive white 
Gaussian noise (AWGN) samples with unit variance at the 
destination and the eavesdropper, respectively. 

III. Optimal Power Allocation 
In this section, we aim to optimize the secrecy performance 
through allocating the relay power Pr, since it affects the 
signal quality at both the destination and the eavesdropper. 
In what follows, we analyze and design the power allocation 
schemes in the sense of maximizing the secrecy outage capac¬ 
ity and minimizing the interception probability, respectively. 


A. Secrecy Outage Capacity Maximization Power Allocation 

Based on the received signal in ([T]), when performing MRC 
decoding at the relay, the channel capacity from the source to 
the relay can be expressed as 

Cs,R = W\og^{l + yR), (5) 

where W is the spectral bandwidth and je = /j||h 5 /{|p 

is the signal-to-noise ratio (SNR). Similarly, according to Q 
and channel capacities from the relay to the destination 
and from the relay to the eavesdropper are given by 

CR,n = ll^log2 (1 + 7d) , (6) 

and 


Cr^E = W log 2 (1 + Je) : 


(7) 


respectively, where je = Prc^r^d 
2 


uH _ 

||h« 


hi?,D 


and 7e = 




hii,. 


PRaR,E 

Then, for the secrecy outage capacity in a DF M-MIMO 
secrecy relaying system, we have the following theorem: 

Theorem 7: Given an outage probability bound by e, 
the secrecy outage capacity is approximated as Csoc = 
W\og2{l+T[aw{Psas,RNE-,PRQ-R,DpNR)) - Wlog^il - 
Pfjafl Ine), when the number of relay antennas is large. 

Proof: Please refer to Appendix I. ■ 

Note that the secrecy outage capacity may be negative from 
a pure mathematical view. Hence, it makes sense to find the 
condition that the secrecy outage capacity is nonnegative. 

For notional simplicity, we let paR^oNR = 
-afl.slne = A ■ ri,Psas,RNR = B, where 
r _ _ ctR.E ine defined as the relative distance-dependent 

‘ paR,DNR ^ 

path loss. Then, the secrecy outage capacity can be rewritten 




Wlog2{l + B)-Wlog2{l + PRAri), B<PrA 
FFlog2(l + PrA) - FFlog2(l + PrAti), B > PrA 


Observing the above secrecy outage capacity, we get the 
following theorem: 
















Theorem 2; If and only if 0 < r; < 1, the secrecy outage 
capacity in such a DF M-MIMO secure relaying system in 
presence of imperfect CSI is nonnegative. 

Proof: Please refer to Appendix II. ■ 

Notice that 0 < r/ < 1 is the precondition for power 
allocation in such a secure relaying system. Given channel 
conditions and outage probability requirements, in order to 
fulfill 0 < r; < 1 , the number of antennas Nr must be bigger 
than — For an M-MIMO relaying system, it is always 

possible to meet the condition of 0 < r/ < 1 by adding more 
antennas, which is one of its main advantages. In what follows, 
we only consider the case of 0 < rj < 1 . 

Based on Theorem 1, the secrecy outage capacity max¬ 
imization power allocation is equivalent to the following 
optimization problem; 

Ji : max Csoc 
Pr 

S.t. Pr < Pmax, (8) 

where Pmax is the transmit power constraint at the relay. 

For the optimal solution of the above optimization problem, 
we have the following theorem: 

Theorem 3: The optimal power at the relay 
is P^ = min Pmax), and the corre¬ 

sponding maximum secrecy outage capacity is 
CZT = W^log 2 (l+min Pmax) - 

IT log 2 (l + min P„,,) An ). 

Proof: Please refer to Appendix III. ■ 

Remark: It is found that when eavesdropper CSI is unavail¬ 
able, it is optimal for the DF secure relaying system to let 
the two hops have the same channel capacity, resulting in 


B. Interception Probability Minimization Power Allocation 


In this subsection, we analyze the optimal power allocation 
from the perspective of minimizing the interception proba¬ 
bility. In general, interception probability is defined as the 
probability of information leakage, namely the probability of 
Cd < Ce- Then, interception probability is equivalent to the 
secrecy outage probability when Csoc = 0 in (fTSl l. Thus, it 
can be computed as 


Po = exp 


2Cd/w _ 2 \ 
PrO!r,e ) ’ 


(9) 


where Cd = W\og 2 (1 + mm{Psas,RNR, PRaR^opNR)) is 
the legitimate channel capacity. 

Then, interception probability minimization power alloca¬ 
tion can be described as the following optimization problem: 


J 2 : min Pq 
Pr 

s.t. Pr < Pmax- (10) 


Since exp(a;) is a monotonously increasing function of x 
and mina;(—/(a;)) is equivalent to maxx{f{x)), J 2 can be 


transformed to the following problem: 

miii{Psas,RNR, PRaR^opNR) 


J 3 : max 
Pr 


s.t. 


Prolr^e 
Pr < Pmax- 


(11) 


By solving the above optimization problem, we have the 
following theorem: 

Theorem 4: From the perspective of minimizing inter¬ 
ception probability, the optimal transmit power at the relay 
P^ belongs to a region 


^ 0 , min ^ 


Psas,. 

paR.D 


P 


and the 


corresponding minimum interception probability is P“'" = 


/ paR,DNR \ 

\ OR.E ) ■ 


exp 

Proof: Please refer to Appendix IV. 

Remarks: It is found that the optimal power minimizing 
the interception probability is not unique. However, from 
the perspective of maximizing the secrecy outage capacity, 
Pr = min(-'^"“®- 


PO:R,D 


P 


namely the upper bound, is 


optimal. Thus, it is better to let P^ = min ^ ’ Fmax 

in the sense of jointly optimizing interception probability and 
secrecy outage capacity. 

C. Asymptotic Characteristic 

In the above, we prove that the optimal relay power P^ 
in the sense of maximizing the secrecy outage capacity and 
minimizing the interception probability is a function of the 
source power Pr. Thus, the source power has a great impact 
on the secrecy performance, as described in Theorem 3 and 4. 
In order to get some clear insights, we carry out an asymptotic 
performance analysis with respect to the source power. 

First, for the secrecy outage capacity in a DF M-MIMO 
secrecy relaying system, there are the following asymptotic 
characteristics: 

Proposition 1: In the low Pr regime, the optimal relay 
power P^ and maximum secrecy outage capacity Cfff^ 
asymptotically approach zero. In the high Pr regime, the 
maximum secrecy outage capacity will be saturated and is in¬ 
dependent of Pr . Furthermore, Cf^ is an increasing function 
of Ps- 

Proof: In the low Pr regime, the optimal relay power 
and the corresponding secrecy outage capacity are reduced 
as PS = and C^T = W^log 2 

respectively. Both of them asymptotically approach zero as Pr 
tends to zero. Otherwise, in the high Pr regime, the optimal 
relay power is limited by Pmax, then Cf^f^ is constant. In 
other words, the secrecy outage capacity becomes saturated. 
In addition, because of 0 < r/ < 1, (7™“ is an increasing 
function of P 5 . ■ 

Second, for the interception probability, we have the fol¬ 
lowing asymptotic characteristics: 

Proposition 2: The minimum interception probability is 
independent of Pr. 

Proof: As proved in Theorem 4, although the optimal 
relay power is a function of the source power, the final 
interception probability is a constant independent of Pr and 
Pr. ■ 
















IV. Simulation Results 


c: 


To examine the effectiveness of the proposed power alloca¬ 
tion schemes for DF M-MIMO secure relaying systems, we 
present several simulation results for the following scenarios. 
We set Nr = 100, W = 10 kHz, p = 0.9 and £ = 
0.01 without extra statements. For convenience, we normalize 
the path loss from the source to the relay as as^R = 1 
and use ur^d and as,E to denote the relative path loss. 
Specifically, ur^e > olr^o means that the eavesdropper is 
closer to the relay than the destination. In addition, we use 
SNRs = lOlog^o Ps, SNRfl = lOlogio^^fl and SNR„,ax = 
10 logj^Q Pmax to represent the SNR in dB at the source, the 
relay and the constraint at the relay, respectively. 


X lo" 



Fig. 1. Comparison of theoretical and simulation results with different outage 
probability requirements. 

First, we verify the accuracy of the theoretical expression in 
Theorem 1 with SNR 5 = SNR/j = 10 dB and ur^d = 1- As 
seen in Fig. [T] the theoretical results are well consistent with 
the simulations in the whole ur^e region with different outage 
probability requirements, which proves the high accuracy of 
the derived results. Given an outage probability bound by 
£, as aR^E increases, the secrecy outage capacity decreases 
gradually. This is because the interception ability of the 
eavesdropper enhances due to the short interception distance. 
In addition, for a given ur^e, the secrecy outage capacity 
improves with the increase of e, since the outage probability 
is an increasing function of the secrecy outage capacity. 

Next, we show the performance gain of the proposed 
optimal power allocation schemes compared to a fixed power 
allocation with SNR 5 = 10 dB, SNRmax = 15 dB and 
aR^E = 5. It is worth pointing out that the fixed scheme uses a 
fixed power Pr = 15 dB regardless of channel conditions and 
system parameters. As seen in Fig.|2] the secrecy outage capac¬ 
ity maximization power allocation scheme performs better than 
the fixed power allocation scheme. Especially in the high ur^d 
regime, the performance of the proposed scheme improves 
sharply, while that of the fixed allocation scheme nearly 
remains unchanged. This is because the legitimate channel 
capacity is limited by the source-relay channel capacity under 
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Fig. 2. Secrecy outage capacity comparison of different schemes. 



Fig. 3. Interception probability comparison of different schemes. 

this condition, but the fixed scheme is regardless of as.fl and 
Ps- In the low aR^o regime, the secrecy outage capacities of 
both schemes approach zero due to r; > 1 , which verifies 
Theorem 2 again. In terms of interception probability, as 
shown in Fig. [2 the proposed scheme also outperforms the 
fixed power allocation scheme. Consistent with the theoretical 
claims, the interception probability approaches zero when 
aR^D is large enough. 

Finally, we check the asymptotic characteristics with 
OiR,D = 1. As shown in Fig. 12 as Pg tends to zero, the max¬ 
imum secrecy capacities with different ur^e approach zero. 
In the large Pg regime, the maximum secrecy outage capacity 
will be saturated, which is in agreement with Proposition 1 
again. From Fig. |5] it is seen that the minimum interception 
probability is independent of Pg. Additionally, the interception 
probability floor becomes higher with the increase of aR^E- 

V. Conclusion 

In this paper, we have first presented a secrecy perfor¬ 
mance analysis for a DF M-MIMO secure relaying system 
with imperfect CSI. We proved that in order to guarantee 
a nonnegative secrecy outage capacity, there is a constraint 



































Vl^log2(l + Psas.fi||hs,fi|p),W"log2 ( 


1 + 


'H r. - ^R,D ^ 


PrOIR^D (v^fl,D + V^l - P^^)- 


||h_R,D 

W^log2(l + Psas,R\\^s,Rf),Wlog^ (l + 


(13) 


PRO!R,Dip\\^R,D\\‘^ + 2 -\/ (1 — p)pTZ{e^hR^D) 
+(l-p)||e^hfi,z.||/||h«.z5f' 


min VFlog 2 (l+ Psas,ij||hs.fl||^), 


Fig. 4. Asymptotic secrecy outage capacity with different olji^e- 



Fig. 5. Asymptotic interception probability with different olr^e- 


on the minimum number of antennas at the relay. Then, by 
maximizing the secrecy outage capacity and minimizing the 
interception probability, we proposed two optimal relay power 
allocation schemes. At last, we revealed the asymptotic char¬ 
acteristics of maximum secrecy outage capacity and minimum 
interception probability with respect to the source power. 


W\og2{l + PRaR^DP\\iiR,D\\^)^ (14) 

Ri W\og2il+inm{Psas,RNR,PRaR^DpNR)). (15) 


where TZ{x) denotes the real part of x. ^ has been replaced 
with ./phji^D + VI - pe in (O. Eq. ([Till follows from the 
fact that p||h_R .dIP scales with the order 0{pNr) as Nr —>■ oo 

while 2v^p(l - p)7^(e^hr^,_D) -f (1 -p)||e^hf_£,|| Vllhfl.cf 
scales as the order 0(1), which is negligible. Eq. (fTSI l holds 
true because of lim = 1 lim = 1. 

namely channel hardening ifT^ . 

Similarly, the eavesdropper channel capacity is given by 


Ce = lElog2(l-l-min(Psas,_R7V_R,PRQ!_R,_E 


1 ^ - 
^R,E I 


iR,D 


Ifl.Cll 


))• 


(16) 

Then, the secrecy outage probability e with respect to a 
secrecy outage capacity Croc can be computed as ([TSb at 
the top of next page, where (ITtI) follows from the fact that 

is distributed with 2 degrees of freedom. 






||hR, 

and (ITFI l holds true since exp 


Ps(^S,rNr ' 
PROiR,E _ ! 

when Nr is sufficient large. Eased on (IlSb 
the Theorem 1. 


approaches zero 
it is easy to get 


Appendix B 
Proof of Theorem 2 


Appendix A 
Proof of Theorem 1 


Based on channel capacities from the source to the relay in 
dUl and from the relay to the destination in (|6]l, the legitimate 
channel capacity can be computed as 


Cd = inin(C5,_R, 

= min ( PPlog 2 (l-f Psas.flUhs.itlP) 


( 12 ) 


W log 2 1 -f PrUr^d 


1 ^ - 
'R,D I 


hfl, 


D 


IR,D\ 


Based on the secrecy outage capacity in Theorem 1, when 


Pr > 


PsOiS,R 

pOiR,D 


we have Caoc = VPlog 2 (l -f PsOts,RNR) - 


lElog 2 (l -f PRpaR^oNRri). To guarantee Csoc > 0, the 
following condition Prc^S.rNr > PRpaR^oNRn must be 
fulfilled, which is equivalent to 0 < r/ < 1 in the case of 


Otherwise, when Pr < the secrecy outage ca¬ 

pacity is changed as Csoc = VPlog 2 (l -f PrPUrpNr) - 
PPlog 2 (1 -f PRpaRpNRTi). Only when 0 < r; < 1, Csoc is 
nonnegative. 

Above all, 0 < r/ < 1 or Nr > — is the 

precondition that the nonnegative secrecy outage capacity 
exists. Therefore, we get the Theorem 2. 
















































Pr{Csoc > Cd — Ce) 


= Pr min i Psas,RNE,PRaR E hfl.-E 

\ llhfl.cl 


= Pr PsO:s,rNr < PrUr^E h|[ E ,, 1 

\\^R,d\\ 


+Pr f PsOiS.R^^R > PrOIR,E hf g < 


( Psas,RNR\ ( PsOiS,RPI 

/ 2 (C'd-C''*°‘=)/^ - 1 \ 


Procr^e 


Appendix C 
PROOF OF Theorem 3 

According to the Theorem 1, when Pr > the 

o in— pafi^D 

secrecy outage capacity Csoc = M^log 2 (l + Psas,RNR) - 
T 4 ^ 1 og 2 (l + PRpaR^oNRri) is maximized when Pr = 
since Csoc is a monotonously decreasing function of 

When Pr < the secrecy outage capacity Csoc = 

fTlog 2 (l + PRpaR^oNR) - fTlog 2 (l + PRpaR^oNRn) is 
an increasing function of Pr under the condition 0 < n < 1 . 
Thus, Pr = is the optimal power at the relay. 

Considering the constraint on the transmit power Pmax 
at the relay, the optimal power at the relay is P^ = 
min f 1 Pmax) ■ Furthermore, by substituting P^ = 

min f 1 P-max j into the expression of secrecy outage 

capacity, we can obtain the maximum secrecy outage capacity 
as shown in Theorem 3. 


First, when Pr < 
equivalent to 


Appendix D 
PROOF OF Theorem 4 

^ PsaS,R 


Cl : max 

Pr 


Pr < min 


, the optimization problem J 3 is 

0!R,DPNr 


PsC(S,R 

POtR,D 


5 Fornax 


Interestingly, it is found that the objective function 
is independent of Pr. Hence, the optimal solution of Gi can 
be an arbitrary value belonging to ( 0 , min ( , Pmax) 1 ■ 


Second, when Pr > 
is reduced as 


Ps<^S,R 


>ptimal so 
^ 0 , min ^ 


Psas.R 


Go '■ max 

Pr 

^ PsCtS,R 

POtR.D 


-, the optimization problem J 3 

Pscis.rNr 

Proir.e 

<Pr< Pro... 


( 20 ) 


^ 2( CD - C , or)/W _ Y 



The optimal solution of G 2 is since the objective 

function is a decreasing function of Pr. 

Thus, the optimal transmit power at the relay is P^ = 

(^0, min Pmax) ■ Hence, we get the Theorem 4. 
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